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Abstract 



' Density-functional theory based global geometry optimization is used to scrutinize the 

O ' possibility of using endohedrally-doped hydrogenated Si clusters as building blocks for con- 



structing highly magnetic materials. In contrast to the known clathrate-type facet-sharing, the 



. clusters exhibit a predisposition to aggregation through double Si-Si bridge bonds. For the pro- 

o 



totypical CrSi2oH2o cluster we show that reducing the degree of hydrogenation may be used to 
control the number of reactive sites to which other cages can be attached, while still preserving 
the structural integrity of the building block itself. This leads to a toolbox of CrSi2oH2o-2M 
monomers with different number of double "docking sites", that allows building network ar- 
chitectures of any morphology. For (CrSi2oHi8)2 dimer and [CrSi2oHi6](CrSi2oHi8)2 trimer 



X 

. structures we illustrate that such aggregates conserve the high spin moments of the dopant 

atoms and are therefore most attractive candidates for cluster-assembled materials with unique 
magnetic properties. The study suggests that the structural completion of the individual endo- 
hedral cages within the doubly-bridge bonded structures and the high thermodynamic stability 
of the obtained aggregates are crucial for potential synthetic polimerization routes via con- 
trolled dehydrogenation. 
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The use of multi-doped endohedral clusters as building blocks for assembly at the nanometer 
scale has attracted much attention due to their promising cage-like geometries and easy-to-tune 
electronic properties.-"- Silicon clusters are particularly suitable for such applications, since en- 
dohedral metal doping was found to produce saturated fuUerene-like cage structures. However, 
the cage geometry in such structures is often stabilized through a strong cage-dopant interaction, 
which leads to quenching of the metal magnetic moment.- Especially interesting for magnetic 
applications is therefore the case of hydrogenated metal-doped Si clusters, where a minimized 
cage-dopant interaction preserves the high magnetic moments of the encapsulated transition metal 
dopants.— 1^ Without doubt, such highly-magnetic building blocks offer great potential for in- 
stance for data storage or spin-based electronic devices.-^^*^ 

Assembling the building blocks into one material structure is, of course, a key problem in 
chemistry and material science. The possibility of building homo- and heterogeneous aggregations 
of MSii6 clusters for different metal dopants M has been theoretically investigated by the groups 
of Balbas^ and Nakajima.^ Robles and Khanna reported that assemblies of CrSii2 clusters may 
have a net spin moment.— However, to the best of our knowledge, no study of potentially highly 
magnetic aggregates of hydrogenated endohedrally doped Si cages has been performed so far. The 
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key question in this is how such saturated clusters can be assembled at all. As we have shown 
in our previous work, Si2oH2o is the smallest Si fuUerene capable of conserving a high spin 
state of an encapsulated 3d transition metal atom, with the H termination simultaneously ensuring 
total saturation of all Si dangling bonds, cf. below. As opposed to this, non-hydrogenated 
MSi2o clusters don't form cage-like structures at all as the Si compensates for its undercoordination 
through strong M-Si interaction. This suggests that a carefully reduced degree of hydrogenation 
(and, therefore, increased amount of unsaturated Si bonds) might be used to control the number of 
bonding sites offered by the cluster. As long as this does not jeopardize the structural integrity of 
the cage-like geometry, such a strategy would lead to a toolbox of monomers with differing number 
of "docking sites", that may offer the possibility to build network architectures of any morphology. 

Scrutinizing this idea through quantitative first-principles calculations is the objective of the 
present study of A/Si2oH2o aggregates. The approach taken is as follows: First, we perform 
density-functional theory (DFT) based global geometry optimizations for the prototypical high- 
spin CrSi2oH2o cage with different number of H vacancies introduced. This indeed validates the 
structural integrity of such "non-ideal" clusters. Second, we also use DFT-based global geometry 
optimization to show that particularly Si-Si doubly bridged dimers and trimers of these clusters 
actually represent most stable and highly magnetic ground- state configurations for a given Si/H 
composition. The energetics thus obtained finally allows us to discuss the thermodynamic feasi- 
bility of such and further polymerization to one-dimensional chains, two-dimensional sheets or, 
eventually, more complicated three-dimensional structures. 

Results and discussion 

Building block monomers 

The starting point of our investigation is the check on the structural integrity of a MSi2oH2o building 
block, when introducing an increasing number of H vacancies to create undercoordinated reactive 
Si sites that would then represent natural "docking" candidates for aggregation. At the minimized 
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3) Si2oH2o b) Si2oH^g 




e) CrSi^^H^g f) CrSi^^H^^ 

Figure 1: Identified ground-state geometries of possible monomer structures, together with their 
spin states and energies relative to their lowest-lying non-cage isomers: a) Si2oH2o, singlet, -3.96 
eV;ii b) Si2oHi9, doublet, -3.18 eV; c) CrSi2oHi9, sextet, -1.23 eV; d) CrSi2oHi8, singlet, -1.95 
eV; e) CrSi2oHi6, singlet, -0.80 eV; f) CrSi2oHi2, singlet, -0.24 eV. 

M-Si interaction in the hydrogenated cage, we do not expect a strong dependence on the actual 
metal M used for doping. This view receives support from test calculations with different metals, 
as well as from the fact that we obtain comparable findings also for the empty Si2oH2o cluster, 
which already by itself adopts a fuUerene-type cage structure, cf. Within the focus on magnetic 
properties, we therefore concentrate in the following on Cr as prototypical dopant atom, which we 
previously reported to yield the highest septet spin state in the CrSi2oH2o cage.— 

Extended configurational sampling shows that one hydrogen vacancy does not change the ge- 
ometrical and electronic structure much, i.e. the identified CrSi2oHi9 global minimum is still a 
(slightly distorted) cage, cf. [Tt, with a rather high sextet spin moment. In order to further an- 
alyze the M-Si interaction and the concomitant magnetic properties, |2] summarizes the total and 
Cr-projected density of states (DOS) of CrSi2oHi9. Consistent with the equally shown 3D spin 
density distribution this confirms that the unpaired electrons are predominantly located on the 
metal dopant, i.e. the latter preserves most of its atomic character as was the case for the ideal 
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Figure 2: Total density of states (DOS) (black solid line) and DOS projected on the metal dopant 
(red dashed line) for CrSi2oHi9 (highest occupied state lies at -5.23 eV, vacuum level is used as 
zero reference). The Cr-specific peak (shown with arrow) indicates that five unpaired 3d electrons 
belong almost exclusively to the Cr atom. The inset shows the spin density distribution within the 
cluster, which resides almost exclusively on the metal dopant. 

CrSi2oH2o cage, 

The next lowest isomer found, which then corresponds to an irregular compact geometry, is 
1.23 eV higher in energy. For the empty Si2oHi9, cf. [ih, this energetic gap to broken-cage ge- 
ometries even amounts to 3.18 eV. In both cases there is thus a clear thermodynamic preference 
for the endohedral cage, which demonstrates that the reduced hydrogenation has almost no effect 
on the structural integrity. What if we thus further increase the amount of H defects? [Hi shows 
the identified ground- state structure for CrSi2oHi8 corresponding to two H vacancies. Again, the 
cage-like geometry is preserved and the energetic gap to the lowest-lying non-cage-like isomer 
is with 1.95 eV quite pronounced. Still, what has changed is the position of the metal-dopant, 
which is now no longer more or less centered within the cage, but located close to the two dehy- 
drogenated Si atoms. This indicates a stronger M-Si interaction of the type found for bare MSi2o 
cages. This perception is confirmed by a spin DOS analysis as the one in|2land also reflected by 
the spin-quenched singlet state of the CrSi2oHi8 cage. 

An intriguing feature that will become central for the aggregation discussed below is the ob- 
tained preference to locate the two H vacancies directly next to each other in the CrSi2oHi8 struc- 
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ture. Energetically lowest-lying alternative arrangements with the two H vacancies e.g. located at 
opposite sides of the cage are by 0.3-0.4 eV less favorable than the paired-vacancy ground-state 
isomer. This structural motif intriguingly prevails if the number of H defects is further increased. 

shows the identified ground-state geometry for the CrSi2oHi6 cluster featuring a total of four 
hydrogen vacancies. Again, the H vacancies are paired and furthermore located at opposite ends 
of the cage. In fact, despite substantial sampling, this was the only H arrangement we could iden- 
tify that yields a stable cage-like geometry at all. All other energetically low-lying isomers found, 
which start above an energetic gap of 0.8 eV to the cage ground-state, correspond to irregular struc- 
tures. Concomitant with the tendency to pair H vacancies, there is thus a tendency to arrange such 
H vacancy pairs as far away from each other as possible. Considering that an agglomeration of 
vacancies at one side of the cage is likely to break it up, the latter tendency is indirectly a conse- 
quence of the energetic stability of the cage, i.e. its resilience to alternative compact geometries. 
Quite nicely, these effectively repulsive interactions between H vacancy pairs then allow to ac- 
commodate an unusually high amount of H defects without jeopardizing the structural integrity of 
the cage. This is also reflected by the obtained ground-state geometry for the CrSi2oHi2 cluster, 
cf. [If, which despite a total of eight H vacancies still maintains an (admittedly distorted, but still) 
endohedral geometry exhibiting four H vacancy pairs. 

Dimerization and trimerization 

The results of the last section reveal the prevalence of the cage-like geometry despite an increasing 
degree of dehydrogenation. One would expect the H vacancy sites formed at the CrSi2oH2o-ii 
cages to then be natural candidates for aggregation, i.e. the docking sites to fuse clusters together. 
With just one H defect per cluster, the obtained CrSi2oHi9 cage only allows for the formation of 
a dimer. However, the resulting dimer shown in|3t is with a computed binding energy of 0.15 eV 
only very weakly bound. Notwithstanding, the geometry shown is the end result of extensive 
sampling runs, i.e. it is the optimum structure found for a given composition of Cr2Si4oH38. In 
particular, this means that the identified dimer is thermodynamically stable against decomposition 
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a) (CrSi,,H^g), 




b) (CrSi^oH^g)^ 



Figure 3: Identified ground-state geometries of cluster dimers: a) (CrSi2oHi9)2, ten unpaired elec- 
trons, 3.83 A intermonomer distance and 0. 15 eV binding energy; b) (CrSi2oHi8)2: twelve unpaired 
electrons, 2.42 A intermonomer distance and 3.46 eV binding energy. 



into any smaller sub-units (under the constraint that in the present microcanonical sampling the 
total number of species in all sub-units together must equal Cr2Si4oH38). Furthermore intriguing 
is the finding that the overall spin moment of the obtained dimer corresponds to ten unpaired 
electrons, i.e. we arrive at a highly magnetic structure. 

As reflected by the 3.83 A distance between the undercoordinated Si atoms in the two cages, 
we attribute the low binding energy particularly to steric constraints that prevent the two clusters 
from further approaching each other without H atoms coming uncomfortably close, cf. [3^. This 
should be much alleviated for the H-vacancy pairs. Indeed, we obtain a significantly increased 
binding energy of 3.46 eV for the dimer formed by two CrSi2oHi8 cages, which after extensive 
sampling again results as ground-state isomer with as large an energetic gap as 6.63 eV to the next 
low-lying structure. As apparent fromlSj), the intermonomer distance is concomitantly significantly 
decreased to 2.42 A and therewith within the normal range found for Si-Si bonds. The resulting 
Si double-bridge type bonding between the two monomers has recently also been suggested by 
Nakajima and coworkers^ as a favorable motif for the aggregation of smaller MSiig clusters. In 
difference to these works, the appealing feature of the present structure is that it affords for a 
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minimized M-Si interaction and yields a high magnetic moment of twelve unpaired electrons. 
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Figure 4: a) Total DOS (black solid line) and DOS projected on metal dopants (red dashed line) 
for the (CrSi2oHi8)2 dimer (highest occupied state lies at -4.67 eV, vacuum level is used as zero 
reference). The inset shows the spin density distribution within the structure; b) Calculated TD- 
DFT optical excitation spectrum of the (CrSi2oHi8)2 dimer ground state, with the insets illustrating 
the dominant Kohn-Sham orbitals behind the characteristic Cr transitions. 



The analysis of the DOS diagram inHJi shows that the overall picture behind this high magneti- 
zation stays the same as discussed for the monomer CrSi2oH2o case before. The distribution of the 
spin density, inset in [4^, further supports this view. Moreover, the optical excitation spectrum in 
|4b, computed within the time-dependent DPT (TD-DFT) linear response formalism-^^^ indicates 
that the optical properties of the (CrSi2oHi8)2 dimer still have a lot in common with those of the 
CrSi2oH2o monomer, and even with an isolated Cr atom.— Similar to the monomer case, several 
well-defined excitations from Cr 4s and 3d orbitals can be clearly discerned, where due to the 
presence of two Cr atoms intensities are just typically higher than in the case of the monomer, and 
some of the peaks are split. 

Encouraged by these findings we have run a global geometry optimization of the linear 
[CrSi2oHi6](CrSi2oHi8)2 trimer, which is at the edge of present-day computational capabilities. 
The resulting ground-state geometry is presented in|5]and indeed corresponds to a highly symmet- 
rical, doubly Si-Si bridge bonded trimer structure with a very high spin moment of 18 unpaired 
electrons altogether, again located predominantly on the three Cr dopants. The energetic gap to 
the next low-lying isomer has increased to an incredible 7.16 eV. We stress again that as a result 
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Figure 5: Ground-state geometry of the [CrSi2oHi6](CrSi2oHi8)2 trimer, featuring 18 unpaired 
electrons. 

of extensive sampling this includes the stability against decomposition into any set of separated 
smaller clusters containing in sum the same number of species. For instance, we have explic- 
itly evaluated by separate global geometry optimization of the smaller fragments that the obtained 
[CrSi2oHi6](CrSi2oHi8)2 trimer is by 8. 15 eV more stable than a set formed of two CrSi2oH2o, one 
CrSig and one Sii2Hi2, or by 9.03 eV more stable than a set formed of two CrSi2oH2o, one Sig and 
one CrSii2Hi2. 

Routes to polymerization 

The dimerization and trimerization results indicate possibility of assembling very stable aggregates 
that largely conserve the magnetic properties of the dopant atoms. Involving cages with higher de- 
gree of dehydrogenation then opens up a pathway to more complicated network architectures, 
e.g. the CrSi2oHi2 cage with four H vacancy pairs might be used for building two-dimensional 
sheets. Another important feature here is that the doubly-bridged Si bonding between the cages 
compensates for the H-vacancy pair induced Si undercoordination in the monomers. Whereas 
CrSi2oH2o-2« monomers with larger numbers of H vacancy pairs were found to keep their struc- 
tural integrity, the increased M-Si interaction still led to a quenching of the dopant spin moment, 
cf. \T\ In contrast, within the dimer and trimer structures the magnetic moment is reestablished 
and contributes to the high total number of unpaired electrons. Indeed, consistent with this un- 
derstanding we have managed to obtain the locally optimized [CrSi2oH8](CrSi2oHi8)6 Si-doubly 
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Figure 6: Locally optimized geometry of a [CrSi2oH8](CrSi2oHi8)6 Si-doubly bridge bonded 3D 
octahedral aggregate, exhibiting 42 unpaired electrons. 

bridge bonded 3D octahedral aggregate shown in |6] with the overall number of unpaired electrons 
reaching an impressive 42. At the implied reminimized M-Si cage interaction within the aggre- 
gates our findings should also not be specific to the Cr dopant used as showcase. In our previous 
work we have established that dopants from the entire 3d transition metal series yield MSi2oH2o 
cages as ground state isomers.-'' Correspondingly, we expect the here suggested dehydrogenation 
route to polymerization to hold for a wider range of in particular magnetic dopants, which would 
then offer a flexible toolbox to engineer electronic or magnetic properties in hetero-aggregates.- 

In our view crucial for a potential synthesis is the structural completion of the individual en- 
dohedral cages in the doubly-bridge bonded aggregates. In this respect, they are real cluster- 
assembled materials and thus differ qualitatively from known Si clathrates which are also com- 
posed of endohedral fuUerene-like building blocks {e.g. MSi2o, ^Si24, MSi28)- Elements of a 
typical Si clathrate usually share facets,—"^ such that all Si valence bonds are saturated by sur- 
rounding Si atoms. They can thus not be separated into structurally complete fuUerene monomers 
and represent rather a space-filling framework than a cluster assembly at the molecular level. In 
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contrast to concomitant clathrate synthesis routes via controlled annealing of bulk material, we 
thus rather expect polymerization strategies via controlled dehydrogenation of formed A/Si2oH2o 
clusters to potentially realize the here suggested aggregates. A central feature for this could be 
their very high stability, which renders e.g. a simple H2 abstraction reaction upon aggregation 
exothermic. For the dimer formation reaction CrSi2oH2o + CrSi2oH2o — > (CrSi2oHi8)2 + 2H2 we 
compute this as 0.35 eV (0. 17 eV per formed H2) released upon aggregation, and similarly 0. 14 eV 
per formed H2 for the abstraction of 4 resp. 12 H2 molecules upon the formation of the trimer and 
3D octahedral aggregate shown in|5]and|6l respectively. This suggests a range of hydrogen chem- 
ical potential to selectively dehydrogenate and build corresponding oligomers, e.g. in solution. 

Conclusions 

In summary, we have sytematically assessed the possibility of using endohedrally-doped hydro- 
genated Si-cages as building blocks for constructing highly magnetic materials. Our unbiased 
first-principles global geometry optimization highlights the structural integrity of the cage geome- 
try of e.g. the CrSi2oH2o fuUerene despite an increasing degree of dehydrogenation. The preferen- 
tially formed H-vacancy pairs then act as natural "docking sites" for polymerization. Depending 
on the amount of H-vacancy pairs of the involved monomers, this yields network architectures of 
any morphology, including linear chains, 2D sheets or 3D structures. 

In contrast to face-sharing Si clathrates, the resulting aggregates represent real cluster assem- 
bled materials, in which structurally complete endohedral Si cages are doubly bridged through 
Si-Si bonds. The latter geometric motif yields a high thermodynamic stability and has also been 
observed in preceding work on smaller MSiig clusters.-^*^ For the present hydrogenated fuUerenes, 
it affords a flexible and controllable molecular assembly, while maintaining the high magnetic 
moments of the dopant atoms inside the constitutent cluster monomers. Taking into account the 
possibility to also stabilize larger cages (such as Si24H24, Si26H26, Si28H28) and combining metal 
dopants having different spin states in hetero-aggregates, this strongly suggests such assemblies 
of hydrogenated endohedral Si-clusters as promising candidates for the construction of highly 
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magnetic nanostructured materials. 



Computational details 

All ground-state total energy calculations in this work have been performed with the all-electron 
full-potential DFT code FHI- aims .-^221 Electronic exchange and correlation was treated within 
the generalized-gradient approximation functional due to Perdew, Burke and Emzerhof (PBE).— 
For comparison the energetic differences between the found isomers were also systematically re- 
computed on the hybrid functional level with the PBEO functional,— without ever obtaining any 
qualitative changes that would conflict with the conclusions deduced from the standard PBE calcu- 
lations. All sampling calculations were done with the "tier2" atom-centered basis set using "tight" 
settings for numerical integrations. The stability of the identified minima has been confirmed by 
vibrational frequency analysis. For the subsequent electronic structure analysis of the optimized 
structures the electron density was recomputed with an enlarged "tier3" basis set.^^ As detailed 
in our previous works systematic convergence tests confirm full convergence of the settings 

with respect to the target quantities (energetic differences between isomers, total electron and spin 
density distributions). To ensure that the obtained geometries indeed represent the ground-state 
structures for all considered systems we relied on basin-hopping (BH) based global geometry op- 
timization,-22i2^ which samples the potential energy surface (PES) through consecutive jumps from 
one local minimum to another. In our implementation 1 1 1 15 ,29 ^j. achieved by random displace- 
ment of atoms in the cluster in a so-called trial move followed by a local geometry optimization. 
A Metroplis-type acceptance rule is used to either accept or reject the jump into the PES minimum 
reached by the trial move. 
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